Quantum-mechanical wavepacket transport in quantum cascade laser structures 
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We present a viewpoint of the transport process in quantum cascade laser structures in which 
spatial transport of charge through the structure is a property of coherent quantum-mechanical 
wavefunctions. In contrast, scattering processes redistribute particles in energy and momentum but 
do not directly cause spatial motion of charge. 



I. INTRODUCTION 

Unipolar quantum cascade (QC) laser devices^ are in- 
traband semiconductor devices in which the transport 
processes and optical (intersubband) transitions which 
give rise to the lasing operation occur only in the conduc- 
tion band of the semiconductor structure. This marks a 
departure from previous interband semiconductor lasers. 
Since the first realization of a QC laser, there has been a 
proliferation of QC laser designsi^ QC laser structures 
are composed of a complicated sequence of semicon- 
ductor layers with different material compositions and 
thicknesses. This sequence is repeated tens or hundreds 
of times giving rise to a periodic structure in the de- 
vice. Considering the complicated layer composition of 
these devices, the nature of the charge transport process 
through these structures is not immediately apparent. 

The original concept by Kazarinov and Suris"^ was 
based on coherent tunneling between neighboring wells. 
It was soon realized however that scattering plays a cru- 
cial role in establishing the nonequilibrium carrier dis- 
tribution, in particular for the depopulation of the lower 
laser level. Thus almost all simulations of transport^i^iSii 
in QC structures assume a semiclassical model, in which 
the transport occurs through scattering transitions be- 
tween energy eigenstates (Wannier-Stark hopping^). In 
such an approach only diagonal elements of the den- 
sity matrix (populations or distribution functions) in the 
Wannier-Stark basis are used, and offdiagonal elements 
are neglected (hence, the term semiclassical). Ref. |^ 
briefly considered a more fully quantum-mechanical ex- 
tension to this semiclassical approach by also including 
offdiagonal elements (coherences) of the density matrix 
in the calculation. This study concluded that quantum- 
mechanical coherences were of limited importance to the 
transport properties. However, coherent effects were ob- 
served experimentally for the electron transfer from the 
injector to the upper laser level in QC structuresi2ii£ 
The importance of coherent effects has also recently been 
stressed by Ref. 11, where a simplified calculation scheme 
is proposed. 

The concept that scattering transitions propel the 
current through heterostructure systems conflicts with 
the standard description of transport in bulk structures, 
where complex Bloch functions carry the current, while 



scattering redistributes the carriers in momentum space 
but does not change their spatial positions. We will 
show that a similar description also holds in quantum 
cascade lasers where the current is carried by quantum- 
mechanical wavepackets, and the scattering only acts lo- 
cally redistributing carriers in energy or momentum. 



II. THE MODEL 

The quantum-mechanical transport theory we use here 
is based on nonequilibrium Green's functions which al- 
low for a consistent combination of scattering and co- 
herent evolution. Quantum-mechanical coherences are 
represented by offdiagonal elements of the G^{E) cor- 
relation function, which is related to the density matrix 
Pak,/3k = ! dE G'^p^{E)/2TTi. The theory is formulated 

with basis states ^Qk(i", 2) = {e^^'^ / ^fA)'4)a{z). Here 
ilia{z) is the envelope function in the growth direction 
z. The wavevector k and the spatial coordinate r are 
two-dimensional vectors in the plane of each semiconduc- 
tor layer (with normalization area A), taking fully into 
account the three-dimensional nature of the structure. 

We divide the total Hamiltonian as H = Ho + ^^scatt ■ 
The free Hamiltonian Ho contains the kinetic energy and 
applied voltage, together with the electron-electron in- 
teraction in a mean field approximation. It is diagonal in 
k, while -ffscatt describes scattering interactions k ^ k'. 
Here we explicitly take into account acoustic phonon 
and longitudinal optical phonon, impurity, and interface 
roughness scattering processes, -ffscatt is treated pertur- 
batively using self-energies in the self-consistent Born ap- 
proximation. For example, for impurity scattering, the 
self-energy has the form 

K^iE) = (C7(q)^i'^(-q))G|,'.k-,(i?) (i) 

/3/3',q 

where y™P represents the impurity scattering 
potentialr^ In the following, we consider both the 
complete form of Eq. |^ with nondiagonal self-energies 
(ND), as well as approximations summarized in Table 
HJ In each case, the system of equations for the Green's 
functions and self-energies is solved self-consistently, 
which is an extensive numerical task^ This results 
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{VccVc^c,) {VcpVpc) {VccVfifi) {VcpV,3'o.') 
DG • • 

ND • • • • 

NDL • - • - 



TABLE I: Scattering potential matrix elements included in 
different self-energy models. • ( — ) indicates scattering poten- 
tial matrix elements included (excluded) from the self-energy. 
The scattering potential V may represent impurity, interface 
roughness, or phonon scattering. The angle brackets represent 
an averaging over the impurity distribution for impurity scat- 
tering, the distribution of thickness fluctuations for interface 
roughness scattering, and phonon modes for phonon scatter- 
ing. DG indicates diagonal self-energies, while ND represents 
nondiagonal self-energies. NDL restricts the scattering matrix 
elements to local terms Vaa for ND. 



the common use of an artificial spatial confinement.) A 
third type of states constitutes the position eigenstate 
basis, using the eigenfunctions of the position operator 
z. Again we start with the Wannier basis, diagonalize 
z, and transform Ho into the new basis. In all bases we 
evaluate the scattering matrix elements directly for the 
respective basis functions. 

Structures: The data presented below was obtained 
using two typical examples: (A) a midinfrared QC laser— 
and (B) a THz QC laser. -'^^ In the calculation we restrict 
the number of basis functions to 9 (5) per period for 
structures A (B). 



in the determination of the full correlation function 
Gp^^^{E), which describes the nonequilibrium state of 
the device. 

Current densities are defined as J = Jo + Jscatt with 



Jo = 



ie{[Ho,z]) 

Vh 



2e M-^ 

Q/3,k 



(2) 

and Jscatt = ieiiHscatt, z]) /Vh. e < is the electron 
charge and V is the sample volume. In Ref. 'T§ it was 
shown that Jscatt provides the hopping current due to 
scattering transitions between the states if one restricts 
to diagonal Green functions and self-energies, i.e. ne- 
glects coherences between the states. Thus this part 
was referred to as scattering current, reflecting the (mis- 
)conception that current flow occurs by a combination of 
both coherent evolution (Jq) and relocation by scatter- 
ing transitions. In the following we will show that Jscatt 
vanishes if coherences are properly taken into account, 
while the entire current is carried by Jo- 

Basis states: The basis functions ipoiiz) can be chosen 
in different ways. In theory, the choice of basis states 
should not affect the physical results. In practice, this 
choice can influence the physical description in several 
ways: i) approximations are always necessary to perform 
the theoretical calculation or to facilitate the numerical 
computation, and the interplay between the chosen basis 
and the approximations may improve or reduce the va- 
lidity of the approximations and hence affect the physical 
result, ii) different choices of basis states can reveal dif- 
ferent physical aspects of the problem. For instance, as 
we show in this paper, working with position eigenstates 
casts light on spatial aspects of the problem. 

In our earlier work,^^ spatially-localized Wannier 
states were used as an orthonormal set of 'ipa{z), which 
can be constructed for the infinitely extended QC struc- 
ture in a straightforward manner and generate a well- 
defined periodic Hamiltonian. A second type of states 
are the Wannier-Stark states, which are the eigenstates 
of Ho- They can be easily obtained by diagonalizing Ho 
in the basis of Wannier states. (This procedure avoids 



III. VANISHING OF J.catt 

In Ref. m we neglected the offdiagonal elements of 
the self-energies, i.e., only terms with a = a' and (3 — 
/3' (diagonal self-energy, DG in Table P) were included 
in Eq. In this restricted theory, we found Jo and 

Jscatt to be similar in magnitude. In contrast, using the 
more rigorous and complete formulation of the theory 
reported here, where all self-energy terms are included 
(nondiagonal self-energies, ND), we find that Jscatt ~ 0. 
This result is basis invariant if the complete nondiagonal 
self-energies (ND) are used, as we checked this explicitly 
(numerically) for all three types of basis sets. Hence the 
total current density is given by Jo- 

The same result has recently been analytically demon- 
strated by Lake and Pandey:^^ In the position eigenstate 
basis, where the matrix Zap is diagonal, the expression 
(16) for Jscatt of Ref. Iisl becomes 



_ 2e y;-^ f dE M-^ / < 
Jscatt - f^y 2^ 2^^""V^a/3,k^, 

k Q/3 



adv 
/3a, k 



Grct v'^ 



x^rct 



Now the term in brackets corresponds to the right-hand 
side of the continuity equation as given in Eq. (57) of 
Ref. IitI which has to vanish. See also Section 2.4 of 
Ref. El 

A second line of argument runs as follows: For mi- 
croscopic scattering potentials which only depend on f 
(but not on momentum p) , the commutator [-ffscatt , z] 
vanishes and hence Jscatt = by its definition. Thus 
Jscatt — holds for a wider class of scattering processes, 
including electron-electron scattering, than studied nu- 
merically here. These analytic arguments demonstrate 
that the previous observation of a finite Jscatt in Ref. ^3 
is an unphysical result caused by the DG-approximation 
in the self-energies for a nonlocal basis of Wannier func- 
tions. 
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FIG. 1: (Color online.) Wavefunctions (modulus-squared) for 
Structure B. (a) Wannier, (b) Position eigenstates. 
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tions, as Vaa represents a transition matrix element be- 
tween states with the same envelope function ipa{z), i.e., 
localized to the same well. This implies that the scat- 
tering processes act only locally in space. They do not 
cause spatial displacement of charge, but act only to re- 
distribute momentum k and energy. 

A key point of our argument was the use of the po- 
sition eigenstate basis which enabled us to demonstrate 
that the dominant scattering processes in Jo act locally 
and therefore do not contribute to the spatial transport 
of charge by Jo- This result is not evident if we work in 
another basis, e.g., Wannier or Wannier-Stark, where the 
wavefunctions are far more delocalized and off-diagonal 
matrix elements Vafs play a significant role. Neverthe- 
less, physical mechanisms, such as the nature of trans- 
port, must be basis invariant. Indeed, as long as the 
full self-energies are taken into account, the expression 
for the currents Jq and Jgcatt are invariant under basis 
transformations. This demonstrates that the same result 
is recovered even in the case of a delocalized basis, where 
the interplay of all types of matrix elements in Table HI 
reproduces the locality of scattering transitions. In con- 
trast, the neglect of certain matrix element combinations 
(such as the DG-approximation) generates spurious non- 
local scattering transition in a delocalized basis. 



FIG. 2: Current-voltage characteristic in position eigenstate 
basis with self-energy models ND and NDL. We display Jo; 
Jscatt ~ in both cases. 



IV. SPATIALLY-LOCAL SCATTERING IN Jo 

Although Jscatt vanishes this does not imply that scat- 
tering processes are absent from the transport process. 
These processes do not appear explicitly in Eq. (|5J but 
they act implicitly in determining G'^{E), and hence 
in driving the current Jq. We now investigate the role 
played by scattering processes in determining Jo- In par- 
ticular, we seek to establish if each scattering event re- 
sults in a spatial displacement of charge, or if a scattering 
process acts locally resulting in energy redistribution but 
with no accompanying spatial transport. 

This issue can be conveniently addressed within the 
basis of position eigenstates. In contrast to the Wannier 
wavefunctions [Fig. Q^a)] , these position eigenfunctions 
are localized to within a single well layer in the struc- 
ture [see Fig. ^h)] and the scattering potential matrix 
elements (a|y|/9) are significantly smaller for a 7^ /9 than 
for a = P- We therefore compare the effect of exclud- 
ing or including certain scattering potential terms in the 
self-energies, see Tabled 

Fig. 121 shows that jf^^ and J^^ are almost identi- 
cal for both samples. This indicates that the potential 
terms (VaaVaa) and {VaotVpf}) constitute the main con- 
tribution to the self-energies in the calculation of J^^, 
our most rigorous formulation of the current. These scat- 
tering terms involve only purely local (in space) transi- 



V. TRANSPORT OF Jo BY COHERENT WAVE 
FUNCTIONS 

As we have just seen, the scattering processes which 
contribute to Jo have a purely local effect which does 
not spatially displace carriers. Hence, we conclude that 
the transport of charge by Jo through the structure must 
arise as a property of the quantum-mechanical wavefunc- 
tion. We can obtain a more concrete visualization of this 
transport mechanism by resolving Jo in energy and space: 

(3) 

The lower panel of Fig. |31 show that the current 
Jo{E, z) flows at different energies in different spatial re- 
gions of the structures. To satisfy the equation of con- 
tinuity, scattering interactions enable energy transitions 
which occur locally as discussed above (an example is de- 
picted by the vertical arrow in the figure). Thus, scatter- 
ing gives rise to a redistribution of electrons in energy but 
does not cause a spatial displacement or spatial transport 
of charge i2i 

Fig. 01 shows the corresponding result for sample A. 
Here the current density is restricted to a narrow energy 
range while crossing the thick barrier (around z ~ —2 
nm) between the injector and the active region, thus 
providing the desired selective feeding of the upper laser 
level. This demonstrates how this representation allows 
for a detailed insight into the microscopic operation of 
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FIG. 3: (Color online.) Lower panel: Spatially and energeti- 
cally resolved current from Eq. ^ for sample B. The wave- 
functions depict some quasi-levels, which are eigenstates of 
Ho at k = renormalized by the scattering interactions. Up- 
per panels: Probability density |(/>nk(2)P (black line) and cur- 



rent (X Re{- 



i'nki{^)/'^i^)} (dashed line) for the complex 



wave functions (f>nii.{z) for k = 0. n corresponds to the largest 
eigenvalue fno{E) for the respective values of E. These ener- 
gies are chosen as local maxima i?max of /„o (E) with FWHM 
AE (~ 1 meV) and the respective occupations have been ap- 
proximated by (7r/2)/„o(-E'max)Ai5. 



the device. 

An especially intuitive form of Eq. © can be ob- 
tained by transforming into an (energy-dependent) basis 
of eigenstates of the hermitian matrix —i/ {271)0^^ k(^) 
which has real eigenvalues fnk{E). While fnk.{E) ex- 
hibits sharp peaks in energy, the eigenstates 4>n]i{z) ex- 
hibit only a weak energy dependence (not shown for 
brevity). Then we find 



Jo{E, ^) = ^ E /"k(i^)Re }^c^l^(z) ^^^^ 



{-ih) d(j)nk{z) 



dz 



(4) 

Thus, the energy-resolved current is carried by wave func- 
tions 4in\i.{z) with occupation densities fn\i{E) (the factor 
2 is due to spin) . This is analogous to the Bloch functions 
and their occupation in standard bulk transport. The 
upper panels of Fig. 13 show that these current carrying 
states 0„k(-z) extend over a range of one period. This 
length scale gives a measure of the coherence length, i.e., 
the distance over which charge is transported by a single 
coherent wave packet. The peaks in density |0nk(-2)P are 
localized in the wells as charge tends to accumulate in the 
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FIG. 4: (Color online.) Spatially and energetically resolved 
current from Eq. Q for sample A. The quasi-levels, which 
are eigenstates of Ho at k = renormalized by the scattering 
interactions are shown for the upper and lower laser state for 
comparison. 



low-potential regions. The currents carried by the wave 
functions (/'„k(-2) are more evenly distributed which re- 
flects the ability of the wave packets to transport charge 
across the structure. Note the nonvanishing divergence of 
current density which reflects the local balance of scat- 
tering processes. Thus, these complex wave functions 
4>n\i{z) can be viewed as generalized states for a nonequi- 
librium system, where the drive by bias and scattering 
processes compensate. 

The currents associated with these wavepackets can 
differ significantly: The state at -29.2 nieV (correspond- 
ing to the upper laser level) has only a small current, but 
a high occupation. In contrast, the state aX E — —41.3 
meV (lower laser level) exhibits a strong current to the 
right, effectively moving electrons away from the lasing 
transition which occurs vertically between these states. 
There are also states exhibiting local current flow to the 
left, e.g., at —22.4 meV. Such considerations may pro- 
vide new tools for the optimization and design of new 
laser structures. 

The complex states 4)n\i{z) are of particular impor- 
tance at level crossings. This is demonstrated in Fig. U 
for a simple superlattice at the resonance between the 
ground state and the first excited state of the neighboring 
well. Fig. |5ja) depicts the Wannier-Stark states, which 
are spread over both wells. The first two eigenvalues 
fno{E) (for k = 0) exhibit a clear peak at the level en- 
ergies [see Fig. EJb)] and the corresponding states, carry 
particles to the right [state 1, Fig. |SJc)] and to the left 
[state 2, Fig.lS^d)]. However, the occupation of the state 
1 is about two orders of magnitude larger thus providing 
a strong current flux over the barrier. 

Our viewpoint of current-carrying coherent wave- 
functions differs essentially from the conception that 
scattering transitions propel the current. This manifests 
itself in two ways: (i) In a scattering transition picture, 
the current would stop immediately once the scattering 
processes stop (TJgcatt drops to zero). However, the oc- 
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FIG. 5: States for the superlattice of Ref. ITgI (a) Dou- 
blet of Wannier-Stark states. (b) Eigenvalues fnk{E) of 
—i/{2-K)G^^j^{E) for k = 0. (c,d) Corresponding eigenfunc- 
tions (pni^iz) at E = — 7meV for the two largest eigenvalues, 
/ik and /2k, respectively. 



swer is that the Wannier-Stark hopping model is derived 
from Eq. ^ by expressing, in a low-order approximation, 
the offdiagonal elements of the density matrix in terms 
of the diagonal elements i^iiSSiS^ Under conditions where 
such an approximation is applicable the Wannier-Stark 
hopping model gives the same results as a full quantum 
transport model (see, e.g., Ref. '53 for superlattices), but 
it should be remembered that the Wannier-Stark hop- 
ping model uses populations which are in actual fact an 
approximation for the coherences. 

Following the concepts of standard quantum optics, 
one is tempted to assume that coherences are created by 
the electric field and destroyed by scattering. However, 
the opposite is true in the quantum treatment of trans- 
port within the basis of Wannier-Stark states, which are 
eigenstates of Ho including the electrostatic potential. 
Here coherences are induced by the scattering processes 
via the non-diagonal self-energies, which are crucial in a 
consistent quantum treatment. 



cupation of current-carrying wave- functions only changes 
on the scale of the scattering time, thus the current flow 
will continue on this time scale, (ii) The current is sen- 
sitive to quantum mechanical phases. E.g., phase conju- 
gation would reverse the direction of Jq, while scattering 
rates would be unaffected. 



VI. RELATION TO WANNIER-STARK 
HOPPING MODEL 

In the Wannier-Stark hopping model, the current is 
carried by scattering transitions between Wannier-Stark 
states and depends only on scattering rates and the pop- 
ulations of these states. The coherences, i.e., offdiagonal 
elements Pak./sk (diagonal in k in contrast to scattering 
transition amplitudes Pak./Jk') of the density matrix, play 
no role in this model and are neglected. There are two pe- 
culiarities in this picture: (i) The Wannier-Stark hopping 
picture of transport through scattering transitions and 
our picture of transport by coherent wavepackets each 
provide the entire current and yet offer such conflicting 
views of the transport mechanism, (ii) If we examine 
Eq. (PI expressed in the Wannier-Stark basis, we see that 
a diagonal density matrix results in zero total current pSi, 
Thus, coherences are central to the transport process, 
even in the Wannier-Stark basis, see also Ref. [23. This 
conflicts with the common Wannier-Stark hopping pic- 
ture. How do we reconcile these contradictions? The an- 



VII. SUMMARY 

We have shown that the entire current through 
QC structures is carried only by complex, quantum- 
mechanical wavefunctions as shown, e.g., in the upper 
panel of Fig. |21 This conclusion follows from two main 
results: (i) The scattering current Jscatt vanishes, which 
directly follows from [iJscatt, ^] = 0, if the scattering po- 
tentials couple to the spatial coordinate of electron states, 
(ii) In the position basis, only the diagonal scattering 
matrix elements Vaa are significant in determining the 
current Jo, and hence the entire current. Such matrix 
elements do not induce spatial translation of the electron 
state within the scattering transition. Thus, the principal 
scattering events which transfer momentum and/or dissi- 
pate energy in the transport process in quantum cascade 
structures occur only through spatially-local transitions 
which do not give rise to propagation of charge. 

The locality of scattering transitions contrasts with the 
simple picture of electrons hopping between energy eigen- 
states due to scattering events, which is a common de- 
scription in heterostructure systems. On the other hand, 
our point of view is analogous to bulk transport, where 
the scattering term in Boltzmann's equation is local in 
space, and the entire current is carried by Bloch states. 
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